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Abstract 

Selective solar absorbers generally have limited effectiveness in unconcentrated sunlight, because of reradiation 
losses over a broad range of wavelengths and angles. However, metamaterials offer the potential to limit radiation 
exchange to a proscribed range of angles and wavelengths, which has the potential to dramatically boost 
performance. After globally optimizing one particular class of such designs, we find thermal transfer efficiencies of 
78% at temperatures over 1,000°C, with overall system energy conversion efficiencies of 37%, exceeding the 
Shockley-Quiesser efficiency limit of 31% for photovoltaic conversion under unconcentrated sunlight. This 
represents a 250% increase in efficiency and 94% decrease in selective emitter area compared to a standard, 
angular-insensitive selective absorber. 
PACS: 42.70.Qs; 81.05.Xj; 78.67.Pt; 42.79.Ek 
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1 Background 

Solar thermophotovoltaic (TPV) systems offer a distinct 
approach for converting sunlight into electricity [1-6]. 
Compared to standard photovoltaics, sunlight is not 
absorbed directly by a photovoltaic material, but is 
instead absorbed by a selective absorber. That selective 
absorber is thermally coupled to a selective emitter, 
which then thermally radiates electromagnetic radiation. 
The key challenge to making such a system efficient is 
achieving a relatively high temperature. Generally, this 
implies high optical concentrations [7]. However, one 
could consider whether there would be another way to 
concentrate heat in the selective absorber-without using 
optical concentrators at all. The key idea here is to 
replace the effect of optical concentration using a differ- 
ent method. 

The most plausible approach to thermal concentration is 
angular selectivity-only allowing light to be absorbed 
within a small range of angles. The reason is that the 
apparent size of the sun is only a small fraction of the 
sky-approximately 1 part in 46,200 [8]. Several researchers 
have considered this in the context of photovoltaics [9] 
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and thermophotovoltaics [6,10]. Metamaterials, such as 
photonic crystals, offer unprecedented control over wave- 
length- and angle- dependent absorptivity. In such systems, 
photon resonances can be tailored to target particular fre- 
quencies and conserved wavevectors to provide pinpoint 
control over thermal emission. Such an approach can be 
applied to create selective solar absorbing surfaces for 
applications such as solar thermal electricity, solar thermo- 
electrics, and solar thermophotovoltaics. The critical figure 
of merit is generally the fraction of incident solar radiation 
capable of being captured as heat. Typically, modest infra- 
red emissivities put strict upper limits on the overall ther- 
mal transfer efficiency possible for the unconcentrated 
AMI. 5 solar spectrum. However, carefully designed 
photonic metamaterials can strongly suppress thermal 
losses in the infrared. 

In this manuscript, we first characterize the performance 
of a standard solar TPV system without angular sensitivity, 
both in the ideal case and with a realistic amount of long- 
wavelength emissivity. We then quantify the improvement 
that can be achieved in a structure with long-wavelength 
emissivity using an optimized angle-sensitive design, as 
illustrated in Figure 1. We subsequently discuss design 
principles for structures with strong angular sensitivity, 
and present calculations on a structure more amenable to 
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fabrication than previous 3D periodic designs [10], consist- 
ing of a 2D array of holes on the surface of tungsten. 

The energy conversion efficiency of a solar TPV sys- 
tem such as in Figure 1 is defined to be [6]: 



rj 



imVin 

ChAc 



(i) 



where 7 m and V m are the current and voltage of the ther- 
mophotovoltaic diode at the maximum power point, C is 
the concentration in suns relative to the solar constant I s 
(usually taken to be 1 kW/m 2 ), and A s is the surface area 
of the selective solar absorber. This system can concep- 
tually be decomposed into two halves: the selective solar 
absorber front end and the selective emitter plus TPV 
diode back end. Each half can be assigned its own effi- 
ciency: r\ t and r\ p , respectively. 

The system efficiency can then be rewritten as: 



n = rj t {T)rip{T) 
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where T is the equilibrium temperature of the selective 
absorber and emitter region. The efficiency of each subsys- 
tem can be further decomposed into its component parts. 
In particular, the selective solar absorber efficiency can be 
represented by [5,11]: 
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where B is the window transmissivity, a is the spec- 
trally averaged absorptivity, s is the spectrally averaged 
emissivity, and a is the Stefan-Boltzmann constant. 

The TPV diode back end efficiency can be represented 
by [6]: 
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Figure 1 Diagram of angle-selective solar thermophotovoltaic 
system. 



where and A E are the effective emissivity and area 
of the selective emitter, respectively. 

2 Results and discussion 

We can begin by considering the situation where absorp- 
tivity for both the selective absorber and emitter is unity 
within a certain frequency range, and S otherwise. The 
ranges for the selective absorbers and emitters are opti- 
mized separately, and the lower end of the selective emit- 
ter frequency range equals the TPV diode bandgap 
frequency co g . If we consider the case of unconcentrated 
sunlight, the limit S — > 0 implies a decoupling between the 
selective absorber and emitter, where the selective absor- 
ber is kept relatively cool to maximize r\ t , while the selec- 
tive emitter acts as if it were much hotter with a bandgap 
frequency co g well over the blackbody peak predicted by 
Wiens law. However, this also leads to declining effective 
emissivity £ E oc 5, and thus A E /A S « IIS. This expectation 
is supported by the numerical calculations in Figure 2 (see 
the Methods sections for details), which demonstrate both 
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Figure 2 For an ideal solar TPV system with unwanted 
emissivity S. a system efficiency versus 5 and b area ratio for 
selective emitter to selective absorber versus 5. 
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that efficiency slowly increases with decreasing 3, while the 
area ratio increases rapidly as 1/3. Clearly the limit where 
3 — > 0 and A E IA S — > oo is unphysical, both because the 
time to establish equilibrium in an arbitrarily large system 
is arbitrarily long, and a perfectly sharp emissivity cutoff 
requires a step function in the imaginary part of the 
dielectric constant of the underlying material However, 
the latter is inconsistent with the Kramers-Kronig rela- 
tions for material dispersion, which derive from causality 
[12]. 

Based on a comprehensive review of selective solar 
absorbers [13], typical spectrally averaged selective solar 
absorber emissivities s are about 0.05 at temperatures of 
approximately 373 K. Assuming 3 = 0.05 as well, this 
implies a maximum system efficiency of 10.5% (T = 720 
K, 7] t = 0.6937, r\ p = 0.1510, A E /A S = 0.75), as illustrated 
in Figure 3a. While a physically relevant result, this effi- 
ciency is unfortunately less than a quarter of the asymp- 
totic efficiency calculated above as 3 — > 0. 

To bridge the gap between performance of solar TPV 
in the cases where 3 = 0.05 and 3 — > 0, we can employ a 
combination of wavelength and angle selectivity. It has 




TPV Bandgap (eV) 



Figure 3 Solar TPV system efficiency: a without angular 
selectivity, b with optimized angular selectivity of functional form 
given in Equation 5. 



been shown in a large number of previous publications 
that absorption can be made to peak at a certain target 
angle or wavevector, over a certain range of wave- 
lengths. While an exact analytical expression is often 
lacking, it generally resembles a top hat function in 
wavelength space, and a local maximum in the angular 
dimension [14,15]. Since local maxima can be approxi- 
mated as inverted parabolas, the analytical expression 
we use is as follows [14,15]: 

s[cd, o) = [i - (e/e m ax) 2 ] [8 + (i - ajn^H], (5) 

where Yl^^^co) is the top hat function, equal to 1 if 
coi <co < co 2 and 0 otherwise. This definition is illu- 
strated in Figure 4 for frequencies within the window of 
the top hat. 

The system efficiency of our angle-selective design was 
determined by inserting Equation 5 into Equation 3, then 
multiplying with the TPV diode back end efficiency of 
Equation 4. Optimizing over the following parameters- 
cutoff frequencies, acceptance angles, TPV bandgap and 
temperature-yields the results in Figure 3b, where the 
maximum efficiency is 37.0% (T = 1, 600 K, r\ t = 0.7872, 
ri p = 0.4697, A E /A S = 0.05). This is 3.5 times higher than 
our previous result, and fairly close to the asymptotic 
limit where 3 — > 0 from before, without the physically 
unreasonable requirement of a perfectly sharp emissivity 
cutoff (which is inconsistent with causality). This result 
also exceeds the Shockley-Quiesser limit for photovoltaic 
energy conversion in unconcentrated sunlight of 31% 
efficiency [8]. Furthermore, as illustrated in Figure 5, 
photovoltaic diodes made from group IV compounds 
such as silicon and germanium have bandgaps that 
would allow for the system to continue to exceed the 
Shockley-Quiesser limit. Finally, the much lower area 
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Figure 4 Schematic diagram of the emissivity as a function of 

angle for all wavelengths. 
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Figure 5 Solar TPV system efficiency as a function of operating 
temperature for germanium and silicon with unconcentrated 
sunlight. Both can exceed the Shockley-Quiesser limit at certain 
operating temperatures. 



ratio A E /A S = 0.05 implies that the angle-selective solar 
absorber illustrated in Figure 1 would serve as a sort of 
thermal concentrator, thus allowing for much less ther- 
mophotovoltaic area to be used compared to previous 
designs in the literature. 

Finally, we consider reasonable metamaterial designs for 
achieving the desired effective emissivity in Equation 5. 
Most structures with nanoscale features on the surface in 
both directions have potential to exhibit strong angular 
sensitivity. The specific structure we examined is a 2D 
array of cylindrical holes in single-crystal tungsten, as dis- 
cussed in [5]. In Figure 6, using numerical techniques 
described in the Methods section, we show that an 
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Figure 6 Emissivity spectra for 2D periodic arrays of cylindrical 
holes in single crystal tungsten at various angles {a = 800 nm, 

r = 380 nm, and d = 3.04 jum. Notice that the average emissivity 
gradually decreases with increasing angle away from normal 
incidence. 



optimal structure with period 800 nm, hole radius 380 
nm, and hole depth 3.04 /mi exhibits decreasing average 
emissivity with increasing angle away from normal inci- 
dence. In particular, at a 75° angle, the average emissivity 
for wavelengths from 400 nm to 2 ftm is 30% lower than 
at normal incidence. Overall, for an absorber in unconcen- 
trated sunlight held at 400 K, the spectrally averaged 
absorptivity a = 0.867, while the spectrally average emis- 
sivity g = 0.073. This results in a projected thermal trans- 
fer efficiency r\ t = 0.750. Such a result compares favorably 
with previously proposed selective absorber designs, such 
as a germanium with a silver back and an anti-reflection 
coating, with a projected thermal transfer efficiency of 
0.678 under identical conditions [5]. Additionally, increas- 
ing the operating temperature to 1,000 K and employing 
100 sun concentration (e.g., with a parabolic trough) yields 
a projected thermal transfer of 0.741; again, above a semi- 
conductor-based design with an anti-reflection coating, 
displaying a thermal transfer efficiency of 0.710 under 
identical conditions [5]. Clearly, suppressing off-angle 
emission with relatively simple structures such as 2D 
arrays of holes in tungsten can give rise to improved spec- 
trally selective performance. Future work should focus on 
modifying these structures to narrow the acceptance 
angles. This approach should allow one to achieve record- 
setting thermal transfer efficiencies for selective solar 
absorbers. 

3 Conclusions 

It was found that although in principle solar thermo- 
photovoltaic systems in unconcentrated sunlight can 
exceed efficiencies of 42%, achieving such performance 
requires suppression of emissivities to unreasonably low 
levels. Conventional materials with undesired emissiv- 
ities of 0.05 display much lower efficiencies of 10.5%. 
However, most of the theoretically allowed performance 
can be restored by introducing angular selectivity of the 
assumed form in Equation 5, with up to 37% overall sys- 
tem efficiency. The system also acts as a thermal con- 
centrator, with receiver areas 20 times larger than the 
emitter areas. Finally, we considered 2D arrays of nanos- 
cale cylindrical holes in single crystal tungsten as a can- 
didate metamaterial for angle-selective operation, and 
found the optimal design parameters to be a period of 
800 nm, a radius of 380 nm, and a depth of 3.04 ftm, 
with a thermal transfer efficiency of 75.0% in unconcen- 
trated sunlight at 400 K. 

4 Methods 

Simulations of electromagnetic properties were con- 
ducted following the same methods as outlined in [5]. 
We employ a finite difference time-domain (FDTD) 
simulation [16] implemented via a freely available soft- 
ware package developed at MIT, known as MEEP [17]. 
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A plane wave is sent from the normal direction and pro- 
pagated through space. On each grid point of a flux 
plane defined at the front and back of the computational 
cell, the electric and magnetic fields are Fourier-trans- 
formed via integration with respect to preset frequencies 
at each time-step. At the end of the simulation, the 
Poynting vector is calculated for each frequency and 
integrated across each plane, which yields the total 
transmitted and reflected power at each frequency [17]. 
The dispersion of tungsten is captured via a Lorentz- 
Drude model [18]. Apart from the approximations of 
material dispersions and grid discretization, these calcu- 
lation methods are exact. 

The emissivity of each structure can be calculated 
from the absorptivity computed above via Kirchhoff s 
law of thermal radiation, which states that the two 
quantities must be equal at every wavelength for a body 
in thermal equilibrium [19]. 

The system efficiency is calculated from numerical 
integration (via the trapezoidal rule) of Eqs. 3 and 4, 
and taking their product as in Equation 2. It can then 
be globally optimized through the application of the 
multi-level single-linkage (MLSL), derivative-based algo- 
rithm using a low-discrepancy sequence (LDS) [20]. 
This algorithm executes a quasi-random (LDS) sequence 
of local searches using constrained optimization by lin- 
ear approximation (COBYLA) [21], with a clustering 
heuristic to avoid multiple local searches for the same 
local minimum. We verified that other global search 
algorithms, such as DIRECT-L [22], yield similar results. 
This ability to directly utilize and compare multiple 
optimization packages on the same problem is provided 
by the NLopt package, written by Prof. Steven G. John- 
son and freely available at http://ab-initio.mit.edu/nlopt. 
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